Stable propagation of an ordered array of cracks during directional drying 
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We study the appearance and evolution of an array of parallel cracks in a thin slab of material that 
is directionally dried, and show that the cracks penetrate the material uniformly if the drying front 
is sufficiently sharp. We also show that cracks have a tendency to become evenly spaced during the 
penetration. The typical distance between cracks is mainly governed by the typical distance of the 
pattern at the surface, and it is not modified during the penetration. Our results agree with recent 
experimental work, and can be extended to three dimensions to describe the properties of columnar 
polygonal patterns observed in some geological formations. 



I. INTRODUCTION 

Cracks appearing when a material shrinks are common 
in everyday life. The most popular examples are proba- 
bly the cracks appearing in paint layers, and those in the 
surface of muds. In the case of paints, there is a superfi- 
cial layer of material that shrinks on top of a substrate to 
which it is attached. For muds, the superficial layer and 
the substrate are the same material, but difference in the 
humidity concentration produces the cracking at the sur- 
face. In these two cases, the shrinking is due to changes 
in the humidity concentration within the sample, but it 
can also be due to the existence of non-uniform tempera- 
ture distributions, which produces stresses and generates 
cracking. The problem of surface fragmentation has been 
studied in the last years, both theoretically and ex- 
perimentally 

There are situations in which cracks appear at the sur- 
face of the material, and penetrate into the sample later 
on. A well known example corresponds to the colum- 
nar fracturing of basaltic lava, seen in many different 
geographical locations ||. A detailed description of this 
problem has only recently been foreshadowed and 
still some points remain obscure. The two dimensional 
equivalent of the columnar structure of basalts corre- 
sponds to the case of a two dimensional material that 
dries (or cools down) starting from a free edge. As in the 
three dimensional case, drying creates internal stresses in 
the material that generate cracks, first appearing at the 
free edge. As the drying front propagates to the interior 
the cracks also do so in turn. There have been different 
experimental realizations of this phenomenon. In one of 
them H the material was a very thin layer of a colloidal 
dispersion placed between two glasses, with one free edge, 
from which humidity can escape. In other experiments 
0, a slurry of A^Oa-water mixture was deposited onto a 
substrate, and a glass was placed a few millimeters above 
it. A low, ultra dry 7V2 gas breeze was injected in the slot 
above the sample. The N2 became saturated with water 
as it passed over the material, and a rather sharp drying 
front propagates with time in the same direction than 
the air flow. A third realization corresponds to the dry- 
ing of thin films of aqueous silica sol-gel, where very nice 
patterns have been observed H. In all these cases the 



propagation of a set of parallel cracks has been observed. 

A detailed theoretical description of this phenomenon 
is lacking, and we will make an attempt in this direc- 
tion for the two dimensional case. In the next section 
we will make a general description of the phenomenon, 
emphasizing the difference with a more standard fracture 
mechanics problem. In Section III we present the calcula- 
tions of elastic and crack energy for the case in which the 
drying (or cooling) front is sharp. In Section IV we dis- 
cuss the conditions under which the crack front (formed 
by the tips of all cracks) is stable, and progresses into the 
material when the drying front advances. In Section V 
we show that the array of cracks generated at the surface 
of the material becomes evenly spaced when penetrates 
the sample, and that the typical width of the stripes is 
determined mainly during the first stage of the process, 
near the surface. In Section VI we briefly comment upon 
the effects of other drying (or cooling) conditions than 
the one studied previously. Section VII contains some 
implications of our work for the study of the three di- 
mensional case. Finally, in Section VIII we summarize 
and conclude. 



II. QUASISTATIC FRACTURE MECHANICS. 
ENERGETIC AND STRESS ANALYSIS 

One remarkable thing about the experiments cited in 
the previous section ^-|| is the fact that a large number 
of cracks penetrate the sample in a coordinated and qua- 
sistatic manner, as the external conditions (the humidity 
profile) change. Then this problem is qualitatively differ- 
ent of a typical problem in fracture mechanics, where the 
propagation of a crack is usually an abrupt phenomenon 
and typically leads to the failure of the material. In stan- 
dard fracture mechanics, stress and energetic analysis are 
two different ways of predicting the evolution of the frac- 
turing process ||. However, the equivalence of both ap- 
proaches is not clear [ ^o| , particularly in cases in which 
the cracks propagate at large speeds, which is almost al- 
ways the case when failure occurs. 

The situation is different in our case. Cracks propa- 
gate only because the external humidity profile changes. 
If the external conditions were stationary, crack advance 
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would be arrested. If we think on the configuration of 
the system as a point in configuration space (the space 
spanned by all coordinates of all particles of the mate- 
rial), at each moment the system is at one minimum of 
the energy landscape. As the humidity profiles changes 
the landscape changes itself, the minimum on which the 
system is located shifts, and the configuration of the sys- 
tem adapts so as to remain in the shifted local minimum. 
This is what we understand as a quasistatic propagation 
of cracks. [p"l| 

Under these conditions the propagation of cracks in a 
two dimensional geometry can be studied by two differ- 
ent but equivalent procedures: Stress analysis consists 
on the calculation of the stress intensity factors |J Kj 
and Kn of the opening and shearing modes at the tips 
of the cracks present in the system. The propagation will 
occur in the direction along which Kn = 0, which coin- 
cides with that for which Ki is maximum. Propagation 
actually occurs only if the energy relieved by the advance 
of the cracks is enough to overcome the fracture energy 
needed to elongate the cracks. In the quasistatic case, 
these two energies will differ only infinitesimally. 

In the energetic procedure to calculate crack advance, 
the total energy (including elastic and crack energy) after 
virtual advances of the cracks is calculated. Cracks will 
actually propagate only if this propagation reduces the 
total energy. Under quasistatic advance the energy re- 
duction during propagation is infinitesimally small, and 
typically, there is only one possible direction of propaga- 
tion for each crack, all other propagation directions would 
produce an increase of the total energy of the system. 

Both approaches are equivalent in the case of qua- 
sistatic advance of the cracks. We will use stress analysis 
or energetic analysis according to convenience in each 
case. 



III. ESTIMATIONS OF ELASTIC AND CRACK 
ENERGY 

We will consider an isotropic and homogeneous ma- 
terial, and assume that linear elasticity can be applied 
|l2"| . It is well known that in this case the material pos- 
sesses only two independent elastic constants fl3[| . Wc 
will take as these two parameters the bulk modulus C 
and the Poisson ratio v. The expansion or contraction 
properties of the material are described by a humidity ex- 
pansion coefficient a, which is formally equivalent to the 
thermal expansion coefficient, i.e., the relative change of 
linear size 6L/L of a piece of material after changing the 
humidity concentration h by some quantity Ah is given 
by 



SL/L = aAh. 



(1) 



The ideal situation we will address consists of a semi- 
infinite two-dimensional sample that is being dried from 



its surface. The drying will be considered to be non- 
homogeneous, and we will model it by a humidity profile 
that is given by some function h(z, t) depending on depth 
z within the material, and time t. The precise form of this 
function will be specified later for different experimental 
situations, but it is important to point out that we do 
not consider the case in which the appearance of cracks 
in the system modifies itself the drying process. Then 
our analysis applies more to cases as that described in 
jjj and possibly ||, but much less to that in || where 
evaporation of humidity through the cracks seems to be 
relevant (see also 14 1). Although in the experimental sit- 
uation the material usually lies on top of a substrate, 
we will study the case in which there is no interaction 
with this substrate. In the experiments of Ref. JtJ this 
is achieved by introducing a layer of some slippery mate- 
rial between the sample and the substrate. Then, all the 
stresses on the sample are originated internally, and are 
due to the existence of a non-uniform humidity concen- 
tration. 

It is known that in an isolated (namely, not clamped) 
piece of homogeneous material placed under a constant 
thermal (or humidity) gradient all stresses vanish in lin- 
ear elastic approximation [ p5[ . Under these conditions 
cracks cannot appear at all, or if already present from 
the beginning its propagation is completely halted. It 
is in fact crucial for the propagation of cracks that the 
humidity gradient it is not constant within the sample. 

We will take the edge of the semi-infinite sample as the 
z = line, and z > in the interior. For convenience, 
we will also refer to the edge as the 'surface', and con- 
sider it to be horizontally placed, in such a way that the 
cracks propagate down the material. At t = the whole 
material is supposed to have a constant humidity concen- 
tration h{). In this situation the material is unstressed. 
Humidity concentration is assumed to decrease with time 
h(z,t > 0) < ho. For any reasonable experimental real- 
ization of the drying process occurring from the z — 
free surface, it is clear that at any time, well inside the 
material we should reach the original humidity concen- 
tration, i.e., h(z — > 00, t) = hg. Then the majority of the 
material is always at constant humidity hg. In this region 
the sample must be unstressed, otherwise it will store an 
infinite amount of energy. Then a boundary condition 
for our problem will be that stresses go to zero as z go 
to infinity. As a simplification of a possible experimental 
situation, we will consider the case of an abrupt drying 
profile (see Fig. |f|), namely h(z,t) = h for z > z (t), 
and h(z,t) = hi for z < zo(t), with hi < ho, and zo(t) 
being an increasing function of time. 

A given set of cracks will always correspond to a local 
minimum of the total energy of the system. The total 
energy is the sum of two well different parts. One is the 
fracture energy, namely, the energy spent in the creation 
of all fractures present in the sample. This is typically 
proportional to the total length of cracks in the material. 
The second part is the elastic energy stored in the sample. 
Let us suppose we have an evenly spaced array of cracks, 



2 



defining stripes of width I, which have penetrated down 
to some distance z, with the drying front being located 
at some position z . 

The clastic energy e c \ stored in the material per unit of 
horizontal length must be proportional to the bulk mod- 
ulus of the material C, and to the second power of the 
typical change in linear density caused by the humidity 
gradient (this change being a Ah, with Ah = (ho — hi)). 
This is at the basis of linear elasticity theory Using the 
fact that I and z — zq are the only two relevant lengths for 
this geometrical configuration [ ft6| ], dimensional analysis 
allows to write down the following expression for e e i 

e el = (aAh) 2 Clg((z-z )/l), (2) 

where g(x) is a dimensionless geometrical function [ jl7| . 
We have determined g numerically and the result is 
shown in Fig. | as a continuous line. We can rationalize 
the general form of the function g(x), considering how 
the elastic energy is distributed along the z direction. 
The total energy e c \ is the integral over z of the density 
of elastic energy Se e \(z). In Fig. |^ we show qualitatively 
the form of Se e \(z) for the cases z <C zq, and z 3> zq. Let 
us consider first the case z < zo (Fig. ||(a)). 5e e \(z) is 
zero for z > zq, since as we already discussed, the ma- 
terial has to be unstressed for z — > oo. In the region 
zq > z > z there is a rather constant energy density, as- 
sociated to the change in humidity concentration, which 
cannot be compensated by a change in density of the ma- 
terial since the material here is attached to the part below 
z p8| . Around the position of the crack front there is 
an increase of the stored elastic energy, which is associ- 
ated to the elastic energy around the tips of the cracks. 
For z < z the elastic energy density goes to zero, since 
here the existence of the cracks has allowed to relieve 
the elastic energy accumulated previously to the cracks 
formation. The linear dependence of g(x) as x — > — oo 
(i.e. for z <C Zq) comes from the energy stored between 
z and zq. Let us consider now the case when z ^> zq 
(Fig. |(b)). The elastic energy becomes independent 
of the crack tips are in a region of material that 

is unstressed. The constant value of g{x) as x — > +00 
comes from the energy stored in the independent stripes 
around the position of the humidity front, which in this 
limit is well behind the crack front. When z ~ zq there 
is a smooth crossover between the two limiting regimes. 

To determine the actual position of the crack front in 
a realistic situation wc will rely on the energetic argu- 
ment. At any time during the drying process, the crack 
front will be located at the position that minimizes the 
total energy of the system. Equation (J2J) gives the elastic 
energy of a set of cracks that have penetrated down to 
the position z. In order to get the total energy, e e ] has to 
be added to the energy cost of creating the cracks. This 
part, when measured per unit of horizontal length, will 
be called the crack energy e c k, and it is simply given in 
terms of the specific energy fracture of the material 77 in 
the form 



e c k = t]zjl (3) 

In order to determine the value of z at which the frac- 
ture front prefers to be located, we have to minimize the 
total energy etot = e e i + e c k- The result we obtain is 
shown in Fig. ^, where we plot the most convenient po- 
sition of the crack front as a function of the parameter 
u = rj/[C(aAh) 2 l] . As we see, for very small value of 
u, x = (z — zq)/1 takes large and positive values, i.e., the 
crack front is located well below the humidity front. This 
is due to the negligible contribution of the crack energy 
compared to the elastic energy in this case. As u increases 
x decreases, crossing zero (namely, the crack front coin- 
ciding with the humidity front) at u ~ 0.66. x tends to 
—00 as u approaches the limiting value uq ~ 1.14. In 
fact, for u > uq the crack energy is so high that cracks 
do not penetrate the sample at all Jlfij ] . 

Under the conditions analyzed in this section, the crack 
front will be located at a position such that (z — zq)/1 
is given by the function plotted in Fig. f|. The crack 
front advances only due to changes in the position of the 
humidity front zq, keeping always z — zq as constant. 

IV. STABILITY OF THE FLAT CRACK FRONT 

The previous analysis has assumed that all cracks pen- 
etrate down to a uniform depth z, and has focused on 
what the value of z is, on energetic grounds. It has to 
be complemented however with a stability analysis of the 
crack front. The necessity of this is clear from the follow- 
ing example. If a material with an array of vertical cracks 
is loaded with a uniform horizontal stress (this situation 
can be thought to be realized in our case if z <C z ), 
there will be typically a single crack that propagates and 
fractures the material. This is a consequence of the fact 
that as soon as a single crack moves forward a small dis- 
tance, the stress at its tip increases, and those at the tips 
of the other cracks decrease. This generates an unsta- 
ble, rapid propagation of a single crack. We will see that 
in our case, this can be compensated by the fact that 
stresses decrease ahead of the humidity front, and this 
can stabilize a flat crack front. 

The same kind of energetic arguments used in the pre- 
vious section will be used to determine the stability con- 
ditions of the crack front. Consider an evenly spaced set 
of cracks, labelled sequentially by an index j, where now 
the tips of the cracks are at vertical positions z 1 , which 
can be slightly displaced from the mean position z, i.e., 
z' J = z + Sij, with Y) , Sj = 0. The horizontal positions 
of the cracks are given by x 3 = jl. The elastic energy of 
this configuration contains a term of the form (^|) , plus a 
correction that can be expanded in powers of Sj/l. The 
first order term in this expansion vanishes, as Y2j $j = 0- 
The second order term can be written in the form 
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3+2 



(4) 



Successive term contain 'interactions' between more dis- 
tant cracks. As the elastic energy is a local quantity, we 
expect \a\ < \b\ < \c\ < ... . Any small displacement of 
the crack front, defined by the quantities Sj can be de- 
composed in a sum of 'normal modes', by going to the 
Fourier representation 



Sj = / dkSkexp(ikx J ), 



(5) 



and the energy decomposes into independent term, in the 
form 

Ae i ~ V] \5k\ 2 {a + bcos(2nk/l) + ccos{^k/l) + ...). 



(6) 

In order for the flat crack front to be stable, Ae e i must 
be positive for any choice of the Si, and thus of the 5k ■ 
This implies that 



(a + 6cos(27rfc/Z) + ccos(47rfc/7) + ...) 



(7) 



must be positive for all k. In the limit of very small k, 
Ae e i is equivalent to a uniform advance of the crack front, 
and then is has to be a positive quantity, as the curvature 
of e e i(z) is always positive (see (||) and Fig. ||). On the 
other hand, It will be proportional to the sum of all coeffi- 
cients in expression (fj), i.e., Ae^j - * ~ a + b + c+ .... Then 
it is clear that if an instability exists for some value of k, 
it will occur at k — n/l, where Ae^ - ~ a — b + c — ... . 
Then we will analyze the stability of the flat crack front 
against a perturbation with k = n/l. We took the equi- 
librium position of a crack front obtained in the previous 
section (plotted in Fig. ||) and calculated numerically 
the quadratic change in energy Ae c i = x(s)(aA/i) 2 C($ 2 /I 
when a perturbation of k = tt/1 and amplitude 5 is intro- 
duced. The dimcnsionlcss function x( x ) can be seen in 
Fig. H (dashed line). We see that x is positive (negative) 
for x greater (lower) than x cr ~ —0.038. In this way we 
obtain that there are two regions in the plot of Fig. ||. 
The one at the left of it cr ~ 1.02 corresponds to a sta- 
ble situation. If the parameters of the system make the 
pattern of cracks to lie in that region, then the time evo- 
lution of the drying process (i.e., the increase with time 
of Zq ) will produce a smooth advance of the crack front, 
keeping always the same value of the distance z — zq to 
the humidity front. At the right of u cr , the pattern is un- 
stable. Should we have one of those patterns at a given 
time, it will immediately propagate forward some of its 
cracks (ideally, one of each two cracks) , in order to reach 
a stable situation. This will imply in particular that some 
cracks will remain halted. The further evolution of the 
crack front will correspond to a new crack pattern with 
less cracks (i.e., with larger I) being propagated. 



V. APPEARANCE AND ORDERING OF 
CRACKS 

Up to now we have assumed that a set of evenly spaced 
cracks exists, and we focused on its stability conditions. 
We will study now how this pattern can appear, starting 
with the process at the surface. Cracks are not expected 
to appear evenly spaced at the surface. In fact, it is 
known in one-dimensional models of surface fragmenta- 
tion |2lJ] (that can be used to represent the first stage of 
cracking of our two-dimensional problem) that the distri- 
bution of fragment length is strongly dependent on the 
presence of small inhomogeneities in the material, mak- 
ing the fragment length distribution broad. But we will 
see that as the superficial cracks penetrate the sample 
they become evenly spaced. We will first discuss how the 
process of nucleating new cracks at the surface is, and 
then argue that the cracks become evenly spaced as they 
penetrate the sample. 



A. The appearance of cracks at the surface 

When the drying front penetrates the material from its 
surface, and before the material gets cracked, the stresses 
can be very simply calculated. The state of the system 
corresponds to the material being completely unstrained 
horizontally (considering the unstrained state as that cor- 
responding to the humidity value ho that occurs well in- 
side the material). Under this condition a uniform hor- 
izontal stress To appears for all z < zq, which is simply 
calculated as To = CaAh. In order for the first crack to 
appear, this value has to overcome the uniform traction 
resistance of the material T r [E2| : 



CaAh > T r . 



(8) 



If this relation is not satisfied, no cracks will appear what- 
soever. Assuming that the relation (||) is satisfied, the 
first crack will nucleate at the surface. This crack pen- 
etrates the sample as long as this penetration reduces 
the total energy of the system. In our case, the crack 
penetrates only down to a distance c?o of the order of 
Zq, where the humidity front is located. The horizontal 
stress at the surface is now zero at the position of the 
crack, and increases as we move away from it, reaching 
the value Fq at large distances of the crack. This means 
that new cracks will nucleate away of the first one, in 
regions where relation (||) is still satisfied. The number 
of cracks nucleated at the surface will be typically the 
minimum number that makes the horizontal stress at ev- 
ery point of the surface to be lower than T r . The typical 
distance between cracks I can then be estimated to be 

I ~ d f(T Q /T r ) (9) 

where the geometrical function f(T /T r ) goes to infinity 
when To — ► T 7 t , and is of order 1 for To T r . Typically, 
I is a few times do- 
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If we stick to the ideal sharp drying front we have been 
studying, the first crack appears for an infinitesimally 
small value of zq, the only restriction is that the rela- 
tion (|h is satisfied. Therefore, since the depth of the 
first crack do is of the order of zq, we should expect a 
very dense set of superficial cracks. However, this un- 
realistic situation is removed if we note that any physi- 
cal drying front will have some typical (finite) humidity 
gradient Vh. At the moment when relation (^) is first 
satisfied (with Ah being now the difference between hu- 
midity concentrations right at the surface, and well in- 
side the material), the first crack will penetrate down to 
a distance 



d ~ Ah/Vh = T r /(CaVh), 



(10) 



and the typical distance between cracks at the surface 
will be a few times this distance. 



B. How the distribution of cracks gets uniform 

The cracks that appear at the surface of the sample 
nucleate at positions that are strongly influenced by the 
presence of small inhomogeneities in the sample . But 
a superficially uneven set of cracks has a tendency to be- 
come evenly spaced as it penetrates the sample. The rea- 
son for this tendency can be understood on an energetic 
basis. For a fixed number of cracks the evenly spaced con- 
figuration corresponds to the minimum of elastic energy. 
It is then clear, using the energetic arguments, that this 
configuration will be approached during the propagation 



process |23 



It is important to clarify the different effects that in- 
homogeneities have on the creation and propagation of 
cracks. For the superficial layer, inhomogeneities are rel- 
evant, and res pon sible for the broad distribution of frag- 
ment lengths pH , since at the beginning all the surface 
is uniformly stressed, and tiny differences in the proper- 
ties of the material will dictate which point of the surface 
will fail first. However, once the array of cracks has been 
defined at the surface, inhomogeneities play a secondary 
role, since new cracks do not appear during penetration, 
and the elastic energy is only very weakly dependent on 
the precise distribution of defects. Then the evolution 
of the pattern is basically the same as if inhomogeneities 
were absent. 

The previous energetic arguments predict the trend to- 
wards evenly spaced cracks as they penetrate the sample. 
Stress analysis allows to re-obtain and make this result 
quantitative. We have calculated the stresses that are 
present in the semi-infinite sample with a sharp drying 
front, in the presence of a particular uneven set of par- 
allel cracks. We took an infinite set of cracks separated 
sequentially by distances l\ and li- 

For this geometrical configuration we have determined 
numerically the direction of maximum opening stress 
(mode I), at the tips of the fractures, and characterized 



it by its angle 9 with respect to the direction of forward 
advance. The first result is that 9 always points in the 
direction of making the pattern more evenly spaced. The 
actual value of 9 is absolutely independent of the values 
of C and r\ of the material, and also on the value Ah of 
the step of the humidity front. It is in fact a quantity 
that only depends on the Poisson ratio v of the material 
and the values of x = (z — zq)/1 and (^i — h)/(h + h)- 
The results for the angle 9, for a material with v = 1/3 
can be summarized as follows: 9 is rather independent of 
x, and within a maximum 10 % error, it can be written 
as 



9(deg)~29(l 1 -h)/(l 1 +l 2 ) 



(11) 



The value of 9 vanishes for l\ = I2, as in this case the 
pattern is actually evenly spaced, and the cracks advance 
straightforwardly by symmetry. The limiting case l\ I2 
gives 9 ~ 29 deg p4| - This is the maximum bending we 
can expect from a given crack, in the process of uniformis- 
ing the widths of the stripes. Note that as a consequence 
of the fact that there are no typical lengths in (pi]), the 
length needed for the pattern to become evenly spaced 
will be proportional to the typical width of the stripes /. 

Once the cracks have started to deviate in order to re- 
duce the elastic energy as much as possible, the further 
detailed prediction of its evolution becomes more diffi- 
cult, since now we should calculate stresses ahead of a 
set of non parallel fractures (or alternatively, the elastic 
energy of this nontrivial configurations of curved cracks). 
However, the main conclusion that the pattern becomes 
evenly spaced in a distance of the order of the stripes 
widths remains lEql. 



VI. DEPENDENCES ON THE 
CHARACTERISTICS OF THE DRYING FRONT 

The sharp drying front is an idealization that is never 
exactly realized in practice. For the experiments of Ref. 
0, for instance, a typical distance over which humidity 
changes is expected. In other cases in which the drying is 
through surface diffusion, the humidity profiles are still 
smoother. As we have already discussed the stability 
reason of a flat crack front in the case of a sharp humidity 
step is due to the rapid reduction of stresses ahead of the 
crack front. Smoother humidity profiles will produce a 
weaker tendency to generate a fiat crack front, and it can 
even occur that a flat crack front is never stable if the 
humidity profile is smooth enough. In this case a pattern 
of bifurcations of the crack front has been predicted [^7) . 
Note that an additional stabilizing factor of the crack 
front exists when the drying is favored by the presence of 
the cracks themselves [jl4| , a case we have not addressed 
here. 
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VII. IMPLICATIONS FOR THE THREE 

DIMENSIONAL CASE 

The same phenomenon we have been discussing ac- 
quires novel characteristics in a three dimensional geom- 
etry. It has a beautiful realization in the geological forma- 
tions named columnar basalts, which can be reproduced 
in a kitchen experiment using corn-starch p8| . When an 
originally hot volcanic lava flow starts to cool down (af- 
ter solidification), the thermal stresses generate cracks 
at the surface, which propagate progressively towards 
the interior of the igneous body. There are many co- 
incidences with the two dimensional case, and also some 
new features ^-^|. Cracks are known to appear at the 
surface and at later times propagate to the interior. At 
the surface the pattern of cracks is rather disordered, but 
becomes progressively ordered as it penetrates the mate- 
rial. In three dimensions, the tendency to order manifest 
in progressive tendency to form a polygonal pattern of 
cracks, in planes parallel to the surface. Qualitatively 
this tendency to order can be understood on the same 
basis as the two dimensional tendency of the 

system to reduce its total energy as much as possible 
during the evolution. It poses however some interesting 
problems, because typically, the expected perfect hexag- 
onal pattern is not reached, but a collection of polygons 
with different number of sides and areas. It has been 
recently demonstrated || that this is a consequence of 
the fact that the patterns starts being disordered at the 
surface, and in the process of minimizing the energy it 
is not able to reach the absolute minimum (namely, the 
honeycomb pattern) but is trapped in a metastable min- 
imum. Based in the arguments for the two dimensional 
case, we can expect also in three dimensions that the typ- 
ical width of the columns is set by the typical distance 
between cracks at the surface, and it is not modified with 
the further penetration. In fact this is what is experimen- 
tally observed, since columns are seen to keep its hori- 
zontal size over distances that reach hundred times the 
width of the columns. For the three dimensional case, the 
cooling through diffusion plus convection (and/or radia- 
tion) at the surface is the most realistic model of cooling 
to consider. Although this type of cooling can lead in 
two dimensions to the instability of the flat crack front, 
this is not necessarily so in three dimensions, since in 
this case all fractures form a connected structure, and 
this generates an additional tendency to stabilize the flat 
crack front. These issues and some additional ones par- 
ticular of the 3D case will be discussed in a forthcoming 
publication. 

VIII. SUMMARY AND CONCLUSIONS 

In this paper we have focused on a two dimensional 
material that cools down or dessicates from one edge. 
We studied the conditions under which a set of cracks 



penetrate the sample during the process. We showed that 
these cracks have a tendency to become evenly spaced, an 
analyze the conditions for which all the tips of the cracks 
form a plane front. Our results provide a quantitative 
framework to analyze recent experiments. Extensions to 
the three dimensional case are expected to provide both a 
qualitatively similar scenario, and also some new features 
that will be fully developed elsewhere. 
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FIG. 1. Sketch of the process studied. There is a sharp 
drying profile at depth zo, with zo increasing slowly with time. 
As this occurs, the crack front (located at z) penetrates more 
deeply into the material. The distance z — zq (which can be 
positive or negative, depending on the parameters) remains 
constant in time. 
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FIG. 2. The functions g(x) (continuous line) and xi x ) (dot- 
ted line) for a material with Poisson ratio v = 1/3. The values 
of x for which x( x ) > are those corresponding to a stable 
flat crack front. 
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FIG. 3. Sketch of the density of energy 8e c \ as a function of 
z for a set of fractures that has penetrated down to z < Zo (a), 
and z > zo (b). The contraction of the stripes is exaggerated 
for clarity. 




FIG. 4. Position of the crack front [x = (z — zo)/l] as a 
function of u = r//[C(aAh) 2 l], obtained by minimizing the 
total energy of the system. The region at the left of u CI cor- 
responds to a stable situation in which all tips of the cracks 
form a stable crack front, whereas at the right of this value 
the flat crack front is unstable. 
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